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ABSTRACT: 

The subject invention discloses a method for non-destructively determining the thickness of 
layers deposited on a substrate by analyzing thermal waves generated in a sample. The methods 
are particularly suited for use with integrated circuit manufacturing. In the subject method, the 
sample is subjected to a focused periodic heat source which generates thermal waves. Either the 
magnitude or phase of the thermal waves generated in the sample are measured. The values 
obtained are normalized relative to a reference sample. The normalized values are analyzed with 
respect to a theoretical model of the sample to calculate the thickness of the unknown layers. In 
an alternate embodiment, thermal characteristics can be determined in a sample as a function of 
depth. The latter approach is useful for nondestructively determining dopant concentrations or 
lattice defects in semiconductor devices as a function of depth beneath the surface. 
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EVALUATING THE THICKNESS OF A LAYER OR DETERMINING 
CHANGE IN THERMAL CHARACTERISTICS WITH DEPTH BY THERMAL 

WAVE DETECTION. 



This invention relates to a method for 
5 determining the thickness of a thin film layer on a 
substrate utilizing a thermal wave detection system. 
In addition, a method is disclosed for developing a 
depth profile of impurities, defects or some other 
depth-varying parameter in a sample. The subject 
10 methods are particularly suited for both detailed 

analysis and production procedures associated with the 
manufacture of integrated circuit devices. 

BACKGROUND OF THE INVENTION 

15 

There is considerable interest in developing non- 
destructive techniques for evaluating subsurface 
conditions* This interest is particularly strong in 
the field of integrated circuit (IQ) manufacturing. 

20 Typically, during the manufacture of an IC paekage, a 
wafer of silicon or other semiconductor material is 
covered with thin film layers. It would be desirable 
to provide a system which is capable of 
nondestructively measuring the thickness of the layers 

25 applied to the semiconductor substrate. 

Another technique used in manufacturing 
semiconductor devices is the diffusion or implantation 
of ions or dopants into the lattice structure of the 
semiconductor. There is a need for a technique for 

30 nondestructively evaluating the concentration levels as 
a function of depth. A suitable depth profiling 
technique could also be used for quantifying lattice 
structure defects, such as vacancies or for measuring 
any other parameter that varies with depth in the 

35 material. As discussed below, the methods of the 
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subject invention satisfy tKe above-stated needs 
utilizing a thermal wave detection system. 

In a thermal wave microscope, thermal features 
beneath the sample surface are detected and imaged by 
sensing the thermal waves that scatter and reflect from 
these features. It is believed that thermoacoustic 
microscopy was first disclosed in U.S." Patent No. 
4,255,971, which is incorporated herein by reference. 
In thermoacoustic microscopy, thermal waves are 
generated by focusing an intensity modulated localized 
heat source at a micoscopic point. As discussed in the 
above cited patent, there are a variety of techniques 
for applying the periodic heat source to the sample, 
for example, an intensity modulated beam of electro- 
magnetic radiation or particle beams. 

Irradiation of a sample with an intensity 
modulated beam of energy results in a periodic heating 
of the sample and in the generation of thermal waves. 
These thermal waves can be measured by a variety of 
techniques depending on vftiich method of detection is 
chosen. One method of detection involves the 
measurement of the oscillating temperature of the 
surface of the sample at the spot of localized heating. 
The oscillating temperature can be measured by placing 
the sample in a photoacoustic cell and measuring the 
pressure oscillations in the cell induced by the 
periodic heat flow from the sample to the gas in the 
cell. (See, "Scanning Photo-Acoustic Microscopy"/ 
y.H. Worig, Scanned Image Microscopy, Academic Press 
London, 1980.) The oscillating surface temperature may 
also be measured with a laser traversing the gas or 
liquid medium in contact with the heated spot on the 
sample surface. This laser beam will undergo periodic 
deflections because of the periodic heat flow from the 
sample to the adjacent medium. (See, "The Mirage 
Effect in Photothermal Imaging", Fournier and Boccara, 
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Scanned Image Microscopy/ Academic Press London, 1980.) 
A third technique for measuring the oscillating surface 
temperature utilizes an infrared detector that measures 
the periodic infrared emission from the heated spot on 
the surface of the sample. (See, "Photothermal 
Radiometry for Spatial Mapping of Spectral and Material 
Properties", Nordal and Kanstad, Scanned Image 
Microscopy, Academic Press London, 1980.) 

Another method for detecting thermal waves 
involves the measurement of the thermal displacement of 
the sample surface at the spot of localized heating. 
Techniques for carrying out the latter method include 
the use of a laser probe or a laser interferometer. 
(See, "Photo Displacement Imaging", Ameri, et al., 
Photoacoustic Spectroscopy Meeting, Technical Digest, 
Paper THA6-2, Optical Society of America, 1981.) 

A third methodology for detecting the thermal 
waves involves the measurement of acoustic signals. 
Acoustic waves are generated by the thermal waves in 
the sample because thermally induced stress-strain 
oscillations are set up in the heated region of the 
sample. These acoustic waves can be detected by a 
variety of techniques including a laser probe (See, 
"Probing Acoustic Surface Perturbations by Coherent 
Light", Whitman and Korpel, Applied Optics, Vol. 8, 
pp. 1567-1580, 1969); a laser interferometer (See, 
"Measurements Using Laser Probes", De La Rue, et al., 
Proc. IEE, Vol. 119, pp. .117-125, 1972); or with an 
acoustic transducer, such as a piezoelectric 
transducer, in acoustic contact with the sample. (See 
U.S. Patent No- 4,255,971, cited above). Any of the 
above-described methods can be used to detect and 
measure thermal waves for performing thermal wave 
imaging and micoscopy. 

In addition to imaging, thermoacoustic microscopy 
can be used for other types of analyses. For example, 
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thermoacoustic microscopy can be used to analyze the 
plate-mode resonant signature of bonded members to 
determine the quality of the bond therebetween. The 
latter technique is disclosed in copending European 
5 patent application no. 83301423. 6, which is 

incorporated herein by reference. As disclosed in the 
present application, thermal wave detection can also be 
used for determining the thickness of a thin film layer 
on a substrate and for obtaining a profile of the 
10 concentration of thermal characteristics in a sample as 
a function of depth. 

Accordingly, it is an object of the subject 
invention to provide a new and improved method for 
determining the thickness of a thin film applied to a 
15 substrate. 

It is another object of the subject invention to 
provide a new and improved method for determining the 
thickness profile of a multilayer thin film structure 
utilizing thermal wave techniques. 
20 It is a further object of the subject invention to 

provide a new and improved method for determining the 
thickness of the topmost layer in a multilayer 
structure using a thermal wave technique - 

It is still another object of the subject 
25 invention to provide a method for profiling as a 

function of depth, a sample that has had its lattice 
structure locally disrupted through the diffusion or 
implantation of foreign ions, such as dopants. 

It is still a further object of the subject 
30 invention to provide a method for profiling, as a 

function of depth, a sample containing imperfections in 
its lattice structure. 

It is still another object of the subject 
invention to provide a method for evaluating a sample 
35 having thermal characteristics that vary as function of 
depth for any reason. 
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It is still a further object of the subject 
invention to provide a new and improved method for 
evaluating the thickness of a layer of material on a 
substrate wherein a thermal wave signal of the sample 
is compared with an expected thermal wave signal 
associated with a reference sample. 

SUMMARY OF THE INVENTION 

In accordance with these and many other objects, 
the subject invention provides for new and improved 
methods for non-destructively analyzing the structure 
of a sample. In order to interpret the results of 
thermal wave detection on a sample, it is necessary to 
construct a mathematical model that provides 
expressions for the temperature at and beneath the 
surface of the sample and for the thermoelastic 
response beneath its surface. If thermoacoustic 
detection is used, the model must also take into effect 
the elastic wave propagation and interference effects 
in the sample. The model must also include all of the 
experimental parameters, such as the character of the 
incoming heat source, the character of the sample, and 
the type of detector used. These parameters must be 
considered since they affect the detected thermal wave 
signals. For example, the magnitude and phase of 
thermal waves generated in material are a function of 
the power and modulation frequency of the beam. In 
addition, the thermal waves are also a function of the 
sample's density, specific heat and thermal 
conductivity. The type of detection system used also 
affects the thermal wave signals. Thus, it is very 
important to consider thermal expansion coefficients 
and the elastic coefficient of the sample wftien thermal 
displacement and thermoacoustic detection is utilized. 
One proposal for developing a model for profiling 
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thermal characteristics in a sample, as a function of 
depth, utilizing a photoacoustic technique, is 
disclosed in applicants book "Photoacoustics and 
Photoacoustic Spectroscopy", A. Rosencwaig, Wiley 

5 Interscience, New York 1980. In order to analyze the 
signals detected as a result of thermal waves, a new 
mathematical model must be developed. Appendix A 
hereinafter presents the derivation of a mathematical 
model of a one-dimensional multilayer system which can 

10 be utilized to carry out the calculations proposed in 
the subject invention, and has been published at pages 
4240 to 4246 in the June 1982 issue of the Journal of 
Applied Physics, Vol. 53. While it is believed that 
the mathematical approach disclosed in Appendix A 

15 represents a fairly accurate characterization of the 
interactions within a sample, it is to be understood 
that the scope of the subject invention is not to be 
limited by the exemplary model. In contrast, in the 
future it is expected that more sophisticated two and 

20 three-dimensional models may be developed which could 
add to the accuracy of the determinations. 

The proposed model provides a basis for 
calculating expected values of thermal wave signals 
based on all the parameters in a given experiment, such 

25 as beam power, the type of detection system used and 
the sample's characteristics. As discussed below, by 
analyzing the experimental data obtained in accordance 
with the subject method, with respect to the model, 
significant information can be developed concerning the 

30 subsurface characteristics of the sample. 

In accordance with one method of the subject 
invention, the thickness of a thin film layer deposited 
on a substrate is determined. In this method, a 
periodic heat source is focused on the uppermost layer 

35 deposited on the substrate. The thermal waves 
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generated are then detected and the value of either 
their magnitude or phase is recorded. 

This value is normalized against a measurement 
taken of reference sample. The reference sample 
consists of a substrate of known construction. The 
normalized values are then compared to values derived 
from a mathematical model corresponding to the 
experimental parameters. In the preferred embodiment, 
the comparison is accomplished using a least-squares 
fitting routine whereby a determination of the 
thickness of the thin film layer can be made. 

The above described technique can be aplied to 
profile the thicknesses of a multilayer deposition on a 
substrate. As disclosed in Appendix A, the multilayer 
. model is expanded recursively to obtain a mathematical 
expression accounting for a plurality of layers having 
different characteristics. Since each layer presents 
another unknown, to solve the equation it is required 
to provide additional experimental data points. 
Accordingly, for multilayer determinations , measurements 
of the thermal wave signals are taken at a plurality 
of modulation frequencies of the heating source. At a 
minimum, the number of test frequencies used must 
exceed the number of unknown layers to be determined. 

In another embodiment of the subject invention, a 
method for generating a depth profile of the 
concentration of dopants, impurities or defects in a 
sample is disclosed. The subject method may be used, 
for example, to determine the concentration, as a 
function of depth, of dopants infused into a structure. 
In the latter method a variation of the multilayer 
mathematical model can be used. More particularly, in 
the mathematical analysis, a portion of the upper 
surface of the sample is divided into hyopthetical 
layers of fixed depth. The unknowns in the equations 
become the thermal conducivity of each layer. In the 


0097473 

- 8 - 

subject method, thermal waves are detected at a 
plurality of modulation frequencies. The number of 
frequencies detected must exceed the hypothetical 
number of layers to be analyzed. The normalized values 

5 are then compared with expected values derived from the 
model • By this arrangement, the thermal conductivity 
characteristics of the sample can be plotted as a 
function of depth. By relating the thermal 
conductivity characteristics to known effects of 

10 impurity levels, the concentration of the dopants can 
be determined - 

In another aspect of the subject invention, a 
method is disclosed wherein the thickness of a layer on 
a substrate can be evaluated. More particularly, in a 

15 production situation, it will be possible to determine 
the expected thermal wave signal associated with a 
desired thickness* Accordingly, by taking proper 
thermal wave measurements, a comparison can be made to 
predetermined values, representative of the desired 

20 sample, for evaluating the thickness of the layer • 
The invention will now be described by way of 
example with reference to the accompanying drawings, in 
which:- 

Pig. 1 is a graphical representation of a two- 
25 layer model for a thin film on a thick substrate. 
Fig. 2 is a graph showing the relationship 
between the phase and magnitude of thermal wave signals 
as a function of depth and frequency, 

Fig. 3 is a graphical representation of a model 
30 for a multilayer system, 

Fig. 4 is a graphical representation of a 
multilayer model of a sample with a nonhomogeneous 
thermal conductivity. 

Fig. 5 is a graph showing the relationship 
35 between the phase and magnitude of thermal wave signals 


- 9 - 


0097473 


as a function of depth and frequency, for a different 
example from that of Fig. 2, and 

Fig. 6 is a graphical representation of a model 
of a sample with a nonhomgeneous thermal conducivity in 
which heating is occurring at a subsurface layer 
x=x 0 =d 0 . 

in a method embodying the invention, for non- 
destructively determining the thickness of layers 
deposited on a substrate by analyzing thermal waves 
generated in a sample, the sample is subjected to a 
focused periodic heat source which generates thermal 
waves. Either the magnitude or phase of the thermal 
waves generated in the sample are measured. The values 
obtained are normalized relative to a reference sample. 
The normalized values are analyzed with respect to a 
theoretical model of the sample to calculate the 
thickness of the unknown layers. In an alternative 
embodiment, thermal characteristics can be determined 
in a sample as a function of depth. The latter 
approach is useful for nondestructively determining 
dopant concentrations or lattice defects in 
semiconductor devices as a function of depth beneath 
the surface. 

Broadly, a method embodying the invention for 
evaluating the thickness of a layer or layers of 
mateial on a substrate by measuring either the phase or 
magnitude parameters of thermal wave signals generated 
by a focused periodic heat source, comprises the 
following steps . 

Focusing the periodic heat source on the uppermost 
layer deposited on said substrate. Measuring the value 
of one of the parameters of the thermal wave signals 
generated in the sample at at least one selected 
modulation frequency of the periodic heat source. 
Comparing the obtained value of the said one parameter 
with a predetermined value of the same one parameter 
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associated with a reference sample, whereby the 
thickness of the layer can be evaluated. The substrate 
maybe a semiconductive material such as silicon. The 
thermal wave signals can be measured with a means that 

5 detects the oscillating temperature of the heated spot 
on the surface of the substrate. The detection means 
may include a photoacoustic cell apparatus, or be 
defined by a laser which senses the periodic heating of 
a medium in contact with the heated spot on the surface 

10 of the substrate. 

Alternatively the detection means may be defined 
by an infrared detector which senses the periodic 
infrared emissions from .the heated spot on the surface 
of the substrate. 

15 Where the thermal wave signals are measured with a 

means that detects the oscillating thermal displacement 
of the surface of the substrate at the heated spot, the 
detection means may include a laser probe, or a laser 
interferometer • 

20 Where the thermal wave signals are measured with a 

means for detecting the acoustic signals that are 
generated by the thermal waves, the detection means may 
include an ultrasonic transducer, or a laser probe, or 
a laser interferometer. 

25 More specifically, to determine the thickness of a 

layer or layers of material deposited on a substrate by 
measuring either the phase or magnitude parameters of 
thermal wave signals generated by a focused periodic 
heat source in a manner embodying the present 

30 invention, the following steps are carried out. 

The periodic heat source is focus sed on the 
uppermost layer deposited on the substrate. The value 
is measured of one of the parameters of the thermal 
wave signals at a plurality of selected modulation 

35 frequencies of the source. The number of modulation 
frequencies selected is greater than the number of 
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layers whose thickness is to be determined. The 
measured value of the parameter of each modulation 
frequency selected is normalized relative to the value 
of the same parameter determined for a reference 

5 sample. The normalized values are compared with 
expected normalized values derived from a model 
depicting the thermal process in the samples, whereby 
the thickness of the said layers can be determined. 
This method of evaluating the thickness of a layer or 

10 layers of material is particularly applicable where the 
substrate is a semiconductive material such as silicon. 
The thermal wave signals are measured with a meas that 
detects the oscillating temperature of the heated spot 
on the spot on the surface of the uppermost layer 

15 deposited on the substrate. The detection means may 
include a photoacoustic cell apparatus, or be 
consituted by a laser which senses the periodic heating 
of a medium , in contact with the heated spot on the 
surface of the uppermost layer deposited on the 

20 substrate, or by an infrared detector vrtiich senses the 
periodic infrared emissions from the heated spot on the 
surface of the uppermost layer deposited on the 
substrate - 

Alternatively the thermal wave signals can be 
25 measured with a means that detects the oscillating 

thermal displacement of the surface of the uppermost 
layer deposited on the substrate at the heated spot, 
such as means including a laser probe, or a laser 
interferometer . 
30 Furthermore, the thermal wave signals may be 

measured with a means which detects the acoustic 
signals that are generated by the thermal waves, such 
as detection means which includes an ultrasonic 
transducer, or a laser probe, or a laser 
35 interferometer. 
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In another embodiment in which the parameter 
measured is the phase of the thermal wave signals, the 
values are normalized by subtracting the value of the 
phase of the reference sample from the phase of the 
multilayer substrate. 

The parameter measured can alternatively be the 
magnitude of the thermal wave signals, and the values 
be normalized by taking a ratio of the value obtained 
from the reference sample and the value obtained from 
the multilayer substrate. 

The comparison of the measured normalized values 
to the normalized values derived from the model is done 
by a least-squares fitting routine. 

To determine the change in thermal characteristics 
as a function of depth in a non-uniform sample having 
impurities or defects therein by measuring either the 
phase or magnitude parameter of thermal wave signals 
generated in said sample by a focused periodic heat 
source in a manner embodying the invention, the 
following steps are carried out. 

The periodic heat source is focussed on the non- 
uniform sample- The value is measured of one of the 
said parameters of the thermal wave signal generated in 
the non-uniform sample at a plurality of selected 
modulation frequencies- The measured value of the said 
parameter at each modulation frequency is normalized 
relative to the value of the said parameter determined 
for a reference sample. The normalized values are 
compared with expected normalized values derived from a 
model depicting the thermal process in the samples, the 
model for the non-uniform sample being characterized as 
divided into a plurality of hypothetical layers each 
having the same thickness, with the number of layers 
being less than the number of modulation frequencies 
selected, whereToy the thermal characteristics of the 
hypothetical layers and the depth variation of these 
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thermal characteristics in the non-uniform sample can 
be determined. The thermal characteristics are the 
thermal conductivities of the hypothetical layers. The 
concentration of impurities or defects in the sample 

5 can be correlated to variations in the thermal 

conductivity such that a profile of the impurities or 
defects can be obtained as a function of depth. 

The measurements of the said one parameter can be 
taken at a plurality of modulation frequencies whereby 

10 the resolution of the said depth variations is 
increased. 

The thermal wave signals can be measured with a 
means that detects the oscillating temperature of the 
heated spot on the surface of said non-uniform sample , 

15 such as detection means which includes a photoacoustic 
cell apparatus, or is constituted by a laser which 
senses the periodic heating of a medium in contact with 
the heated spot on the surface of the non-uniform 
sample, or by an infrared detector which senses the 

20 periodic infrared emissions from the heated spot on the 
surface of said non-uniform sample. Furthermore, the 
thermal wave signals can be measured with a means that 
detects the oscillating thermal displacement of the 
surface of said non-uniform sample at the heated spot, 

25 in which case the detection means may include a laser 
probe, or a laser interferometer. 

Alternatively, where the thermal wave signals are 
measured with a means which detects the acoustic 
signals that are generated by the thermal waves, the 

30 detection means may include an ultrasonic transducer, 
or again a laser probe or a laser interferometer; 

If the parameter measured is the phase of the 
thermal wave signals, the values are normalized by 
subtracting the value of the phase of the reference 

35 sample from the phase of the non-uniform sample. 
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If the parameter measured is the magnitude of the 
thermal wave signals, the values are normalized by 
taking a ratio of the value obtained from the reference 
sample and the value obtained from the non-uniform 
sample. 

The comparison of the measured normalized values 
to the normalized values derived from the model can be 
done by a least squares fitting routine. 

In order to analyze the thermal wave signals 
detected when a sample is tested, a mathematical model 
must be developed. Referring to Fig. 1, a two-layer 
representation for a thin film layer 20 deposited on a 
thick substrate 22 is illustrated. Using such a 
representation, mathematical expressions can be derived 
far explaining the effects of thermal waves. 

As discussed above and in the previously cited 
references, thermal waves can be detected by measuring 
variations in temperature or thermal displacement at 
the surface or by detecting the thermoacoustic signal 
generated by the thermal waves. In any case/ the 
detected signal will be periodic in response to the 
amplitude nodulated incoming heat source. These 
signals are quantified by correlating either the 
magnitude of the output signal or its phase in relation 
to the incoming modulation signal. 

The output signal detected will depend on a 
variety of parameters, including the frequency of the 
modulation signal the method of detection and the type 
of material being tested. Because there are so many 
factors affecting measurement, an absolute standard 
cannot be conveniently used. Accordingly, measurements 
of the thermal waves must be normalized against a 
reference sample. As discussed below, the reference 
sample is defined based upon the, type of information 
which is being sought. 
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As discussed above, the normalized values of 
either the phase or magnitude of the thermal wave 
signals are used to determine the desired information. 
This result is possible because the magnitude and phase 

5 of the thermal wave signal varies with the modulation 
frequency for a particular ratio of two thermal 
conditions. The graphs displayed in Fig. 2 plot 
normalized values for the magnitude and phase 
parameters, in a two-layer model illustrated in Fig. 1 

10 where the thermal conductivities K are related as 
K 2/ K l" 4 * In *-he graphs, the X axis is frequency 
normalized in units of d/p where d is the thiclcness of 
the top layer and ^i is the thermal diffusion length 
given by 

15 p = (— ) 1/2 

where J> is the density, C is the specific heat and CO 
is the beam modulation frequency in radians per second. 

20 The curves illustrate how the normalized parameters of 
magnitude and phase vary over both the depth of the 
layer and frequency of the modulation beam. 

Turnign now to the method of the subject 
invention, a value for either the phase or magnitude of 

25 the thermal waves, needed to achieve normalization, is 
obtained by focusing a periodic heat source on a 
reference sample. In the situation where the thiclcness 
of a single layer deposited on a substrate is to be 
determined, the reference sample can be the uncoated 

30 substrate . Since the phase and magnitude of the 
. detected signals are not dependent upon the thiclcness 
■ of a uniform reference sample, any sample corresponding 
to the substrate to be tested may be used as a 
reference sample. 

35 A periodic heat source 24 is then impinged on the 

upper surface of thin layer 20. As discussed in 
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Appendix A, the thermal wave generated will be 
scattered and reflected from the boundary between the 
thin layer 20 and the substrate 22. This effect can be 
observed by noting the variation either in the 

5 magnitude or the phase of the thermal waves detected* 
The magnitude or phase values measured are then 
normalized, to derive a value which is representative 
of the layer to be analyzed* As is well-known, 
magnitude values can be normalized by calculating the 

10 ratio between the reference value and the sample value. 
Since phase measurements are periodic, the 
normalization of these values can be obtained by 
subtracting the refernece phase from the test sample 
phase* A measurement of either the magnitude or the 

15 * phase can be used in the method of the subject 
invention. In many cases, the measurement of 
normalized phase will be preferred since phase 
measurements are dependent on fewer experimental 
variables than are the magnitude measurements. 

20 Once the paramter values have been normalized, 

they are compared with expected -normalized values 
derived from the model which is being utilized. One 
suitable model is disclosed in Appendix A. The 
mathematical model takes into account all experiment 

25 parameters, and is arranged in an expression which 
varies as a function of depth. 

In order to obtain an unambiguous determination of 
the thickness H d tt of layer 20, it is preferable to make 
measurements on the sample at a plurality of beam 

30" modulation frequencies. Thereafter, the experimental 
data obtained is applied to a mathematical equation 
defining the model and the equation is solved for "d*« 
Preferably, the equation is solved using well-known 
least-squares fitting routines Which optimize the fit 

35 between the experimental results and the theoretical 

model* The accuracy of the results will be enhanced by 
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providing additional experimental data points derived 
by obtaining a plurality of information at various beam 
modulation frequencies . 

The above described technique provides for a 

5 method for yielding a quantitative measurement of the 
thickness of an uppermost layer deposited on a 
substrate. In some manufacturing situations, it is not 
necessary to caluclate the actual thickness of the 
layer deposited but rather to determine if that 

10 thickness corresponds to a desired thickness. Stated 
differently, in a manufacturing setting, it may be 
necessary to determine only if the thickness of a 
particular coating corresponds to the desired thickness 
and is therefore satisfactory, rather than having to 

15 calculate the actual thickness of the coating. 

In the latter situation, it may be unnecessary to 
compare the data measured during testing with 
theoretical data derived from a normalized model. 
Rather, the experimental results can be compared to a 
. 20 predetermined expected result, corresponding to a 
sample having a layer deposited with the desired 
thickness. Thus, the subject invention contemplates a 
method wherein the experimentally derived normalized 
parameter is compared to a predetermined normalized 

25 parameter for evaluating whether the thickness of the 
layer in question corresponds to the predetermined 
thickness layer. As with the first described method, 
in order to eliminate any ambiquity, it would be 
desirable to take a plurality of measurements at 

30 various modulation frequencies. 

The latter method may be particularly suited for 
integrated circuit manufacturing technique wherein 
layers of material are depositing on a silicon or other 
seraiconductive substrate. After arriving at a 

35 predetermined value corresponding to a satisfactory 

construction, each manufactured IC can be scanned and 
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measured. If the measured values correspond to the 
predetermined values, the coating thickness can be 
classified as satisfactory. 

In order to obtain a thickness profile of a multi- 
layer thin film structure, a multilayer model must be 
developed. More particularly, and referring to Fig. 3, 
the number of layers in the model must be equal to the 
number of layers in the structure to be evaluated. 
Using a reitterative version of the two-layer model 
discussed above, an expression representing the 
thicknesses of the various layers can be calculated. 

In this analysis, the thermal parameters (Ki-Kn)- 
of each layer is known. In order to calculate all of 
the unknown thicknesses (d^d^, it is necessary to 
experimentally provide a number of data points vftiich 
exceeds the number of unknown thicknesses. 
Accordingly, in the method of the subject invention, 
the values of parameters are determined at a plurality 
of selected modulation frequencies herein the number 
of modulation frequencies selected is greater than the 
number of unknown layers. In this situation, wherein 
the thickness of a plurality of layers is to be 
determined, the normalization of the values obtained 
will be relative to the underlying uncoated substrate. 
As in the previous example, this reference sample can 
be either the actual sample prior to deposition of the 
layers or another sample having a similar structure and 
thermal characteristics. 

The depth of penetration of the thermal waves for 
effective imaging is dependent upon their wavelength. 
More particularly, longer wavelengths, and hence lower 
modulation frequencies will provide images to a greater 
depth. Thus, in multilayer analyses, the frequencies 
selected must include wavelengths sufficient to provide 
information as to the lowest layer of interest. 
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In order to determine the thickness of only the 
topmost layer, in a multilayer structure, where all the 
other layer thicknesses are known, multiple 
measurements, as required in the previously described 

5 method are not required. More particularly, when the 
other thickness layers are known, their values can be 
applied to the theoretical multilayer model with only 
the thickness of the uppermost layer being unknown. In 
this instance, the reference sample used for 

10 normalization will be the multilayer structure prior to 
the topmost layer being deposited. By this arrangement, 
the experimental values derived will represent only 
that of the uppermost layer. As in the determination 
of a single layer thickness, measurements at additional 

15 modulation frequencies will reduce the likelihood of 
ambiguities in measurement. 

By using the multilayer model in a different 
approach, it is possible to profile a sampled thermal 
characteristics as a function of depth. This technique 

20 is of extreme interest in .the field of integrated 

circuit manufacture where it is desirable to examine 
the lattice structure of semiconductor materials, such 
as silicon. More particularly, manufacturers, 
frequently diffuse or implant ions or dopant materials 

25 in the silicon structure in order to effect its 

semiconductive properties. The subject method can be 
utilized to develop a profile of the concentration of 
these dopants as a function of depth from the surface 
of the material. Another use for the subject method 

30 would be to analyze defects, such as dislocations or 
vacancies, in the lattice structure as a function of 
depth. 

As can be appreciated, in the above discussed 
multilayer model, determination of layer thicknesses 
35 was made based upon prior knowledge fo the thermal 
characteristics of each material layer deposited. 
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Thus, while the thicknesses were unknown, all other 
thermal characteristics were known. In this embodiment 
of the subject method, the sample is treated as a 
mulilayer structure having unknown thermal 
characteristics. Further, the model is divided into a 
plurality of hypothetical layers of a known thickness. 
Thus, the unknowns in the mathematical model become 
thermal characteristics such as the thermal 
conductivity, while the thickness of the layers are 
known. As discussed below, by calculating the thermal 
characteristics as a function of depth, a depth profile 
of dopant concentrations can be derived. 

In accordance with the subject method, a periodic 
heat source is focused on the nonuniform sample. 
Measurements of one of either the magnitude or phase 
parameters of the thermal wave signal are taken at 
selected modulation frequencies. Similar to the 
multilayer method, the number of frequencies used must 
exceed the number of hypothetical layers being 
investigated. As can be appreciated, resolution of the 
depth profile can be increased by increasing the number 
of hypothetical layers. However, since an increase in 
the number of hypothetical layers requires an increased 
number of test frequencies, the length of time for an 
examination will increase. Thus, in a manufacturing 
situation, the need to maximize the resolution of the 
depth profile can be balanced against the time 
necessary to carry out the measurements. 

As in the previous methods, the values obtained 
must be normalized against a reference sample. In the 
instant method, the reference sample is characterized 
by a uniform or nontreated material. By this 
arrangement, the effects of the dopant or lattice 
irregularities can be directly analyzed. 

Similar to the calculations described for the 
multilayered determination, the measured data is 
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analyzed with a least-squares fitting routine relative 
to the multilayer model. As pointed out above, in this 
analysis, the thicknesses of the layers are known, 
while the thermal characteristics of the layers are 

5 unknown. Once the thermal conducitivites of the layers 
are calculated, it is necessary to obtain a correlation 
between the thermal conductivity of the lattice with 
respect to the concentration of dopants or defects in 
the lattice. Stated differently, a relationship must 

10 be established between the value of the thermal 

conductivity and the level of dopants or defects in 
the lattice. Such a correlation may be obtained from 
calibration experiments or from direct measurements of 
thermal conductivity for different concentrations of 

15 dopants or defects. 

Thus, by this method, values can be derived which 
represent the concentration levels of dopants or 
impurities in each hypothetical layer. This 
information provides a profile, as a function of depth 

20 of the desired characteristic. 

In summary, there has been disclosed new and 
improved methods for determining the thickness of a 
layer of material deposited on a substrate by studying 
the phase or magnitude parameter of the detected 

25 thermal wave signals. In the subject method, a 

periodic heat source is focused on the sample. A 
measurement is made of one of the parameters of the 
thermal wave signal. The measurements are taken at a 
plurality of beam modulation frequencies with the 

30 number of frequencies selected being greater than the 
number of layers whose thicknesses are to be 
determined. The values are normalized relative to the 
values obtained from a reference sample. The 
normalized values are then alayzed with respect to a 

35 mathematical model of the multilayer system. The model 
represents a set of equations taking into account all 
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phenomena associated with the thermal wave generation 
within the sample. Using a least square fitting 
routine, the unknown thickness of each layer can be 
calculated. In an alternative method of the subject 

5 invention, a profile of the thermal characteristics of 
a sample caused, for example, by dopants or lattice 
defects can be determined as a function of depth. In 
the latter case, the theoretical model is divided into 
a number of hypothetical layers of known thickness. 

10 Using a least-squares fitting routine on data obtained 
at various modulation frequencies, the thermal 
characteristics, and in particular the thermal 
conductivity of the hypothetical layers can be 
determined. The thermal conductivities can be related 

IS to the concentration of dopants or lattice defects to 
provide an accurate depth profile of the sample. 
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APPENDIX A 
THERMAL-WAVE DEPTH-PROFILING: THEORY 

SUMMARY 

A one -dimensional model for thermal-wave 
depth-profiling provides expressions for the temperature 
at the surface of the sample and for the thermoelastic 
response beneath the surface. The model shows that 

5 elastic wave interference effects produce significant 
differences between samples with mechanically free and 
constrained surfaces, and that thermal-wave images of 
thermal conductivity variations are obtainable from the 
thermoelastic signal only if the front surface is 

10 mechanically free. The case of subsurface heating 
shows that for heating occurring at depths of more 
than a few thermal diffusion lengths, the thermoelastic 
signal becomes independent of thermal conductivity 
variations- This has important implications for thermal- 

15 wave image range and resolution. 
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It has been suggested that one can, with the 
photoacoustic technique/ obtain information about the 
optical, and /or thermal characteristics of a sample 
as a function of depth beneath its surface. 1 This 
nondestructive depth-profiling is a unique capability 
of photoacou sties because of the critical damping of 
the thermal waves that are generated in the process. 
Although there has been some experimentation on 

1 2 

photoacoustic, or thermal -wave, depth -prof iling ' , 

this capability has not been extensively exploited, 

primarily because of the lack of adequate theoretical 

models. There have been some .theoretical analyses of 

3 4 

thermal-wave depth-profiling ' but these appear to 
suffer from several inadequacies. The advent of 
thermal-wave imaging and microscopy, with its 
capabilities for subsurface flaw detection and for 
three-dimensional imaging of semiconductor dopant 
regions, has generated a need for a rigorous, yet 
usable, analysis of thermal-wave depth-profiling. 
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Thermal-wave imaging of subsurface features has been 
performed with gas-microphone photoacous t ics 5 ' 6 f with 
photothermal techniques 7 ' 8 and with piezoelectric 
photoacoustics 9 . These initial experiments were all 
performed at low modulation frequencies, and since the 
resolution in a thermal-wave image is set both by the spot 
size of the illuminating beam, and by the thermal 
wavelength 10 ' 11 truly microscopic* thermal-wave imaging is 
possible only at high modulation frequencies. Such high 
resolution thermal-wave images are presently possible only 
with piezoelectric photoacousVic detection, and have been 
demonstrated both with laser beam 12 and with electron 
beam scanning 13 * 14 . 

Experimentally, thermal-wave imag i ng per formed wi th 
piezoelectric detection bears an apparent resemblance to 
thcrmoelast ic ultrasonic imaging. White 15 first 
investigated the phenomenon of thcrmoelast ical ly generated 
ultrasonic waves, and since then there have been several 
applications of the phenomenon, both for materials 
evaluation 16 ' 17 and for ultrasonic imaging 18 ^ 21 . 
Thermal-wave imaging is, however, not the same as 
thermoelastic ultrasonic imaging. In thermal-wave imaging 
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il is the interactions of the thermal -waves with the 
thermal features in the sample that produce the main image 
features, while in thermoelas t i c ultrasonic imaging it is 
the interactions of the acoustic waves with the elastic 
features in the sample that produce the principal image 
features. Thus, as defined above, a thermal-wave image can 
be obtained either through a measurement of surface 
temperature variations, as in the cases of gas-microphone 
photoacoustics and photothermal detection, or through a 
measurement of the subsurface thermoelast ical ly generated 
ultrasonic waves. The comnon element between these three 
quite different detection techniques is that the 
thermal -wave image arises from the interaction of the 
thermal waves with the sampled thermal features. Thus 
such aspects as depth-profiling and the imaging of dopant 
regions in semiconductors are unique Matures of 
thermal-wave imaging and are not possible in thermoelastic 
ultrasonic imaging. The depth-profiling capability arises 
as we have already stated, from the critically-damped 
nature of the thermal waves, and dopant regions in 
semiconductors are visible to thermal waves because they 
undergo noticable changes in thermal conductivity while 
still having essentially unaltered elastic properties. 
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The generation of thermoelastic waves has been treated 
theoretically by several authors 15 ' 20 ' 23 as has the case 
of beam-induced thermal transients 24,25 . However, 
neither the thermal-wave imaging aspects in general, nor 

5 the depth-profiling aspect in particular, have been fully 
addressed in a model that considers both the effects on 
surface temperature and on thermoelastic response. The 
theory presented here attempts to present such a model, one 

lo that is applicable to gas-microphone, photothermal and 

piezoelectric thermal-wave experiments. In the latter case 
the important role of sample elastic properties is 
presented in the Initial parts of the theory. However, for 
the sake of a more simplified analysis and a cleaner 

15 Physical interpretation of the results, the theory is fully 
developed for the special case where the sample has uniform 
elastic properties. By so doing, we more clearly show the 
inherent similarities and differences between thercial-wave 
imaging performed by measuring surface temperatures, as in 
2o gas-microphone or photothermal experiments, and 

thermal-wave imaging performed by measuring thermoelastic 
response in the sample. The special case of a sample with 
nonuniform thermal but with uniform elastic properties is 
moreover appropriate for what is probably the most 
25 important application of thermal-wave depth-profiling, 
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that is, the profiling of dopant concentrations in 
semiconductors. In addition, the theory can be used to model 
the potentially powerful technique where both the phase as well 
as the frequency are varied. We have also considered the case 
5 where the energy absorption and subsequent heeling occurs al 
some depth beneath the surface of the sample. The results 
obtained from the analysis have important implications to 
thermal-wave resolution and image range. 
II. Theory 

1° Consider a plane heat source with sinusoidal time dependence, 
Q Q e icjt , at the surface, x=x o =0, of a semi-infinite elastic 
body. At any practical frequency, we will be dealing with the 
uncoupled thermo-elast ic problem and for this we need to first 
solve the homogeneous one-dimensional heet equation, 

V 

for the temperature, T, as a function of position, x, subject to the 
boundary condition, 

.^11 = Q o (2) 

dx Jx=0 ° 

In Eq. (l) f q is a thermal wave vector defined by 

q " (1 '^") (3) 
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where p is the density, C the specific heat, ond k the thermal 
conductivity and ji is commonly referred to as the thermal diffusion 
1 ength. 

Let us consider the sample to be a thermally inhomogeneous system 
5 where q depends on x. The resulting solution for T for this case is 
in general quite complicated. A straightforward and illuminating 
approach to this problem is to model the inhomogeneous portion of the 
material with a system of N plane homogeneous layers, as depicted in 
Fig. 4, each of whose thermal properties correspond to some average 
, lo over layer thickness of the material's thermal properties- The 
temperature within the nth layer is then given by 


Tfx) = A n e" < 'n x + Be q n x , (4) 

n n n . 


■where 


15 and the solutions within adjacent layers are related through the 
boundary conditions, 

T n|x=x n = T n+l|x=x n 

(5) 

, 1*\ - dT nti| 

*n 3x Jx=x n *n+l dx~ l x=x n 

let us define a characteristic thermal-wave impedance of the nth layer 
20 Zn = *n<Jn < 6 > 

and a thermal-wave input impedance, Z^ n f to the nth layer by 

,in . *n dx I /,\ 

Z n T~- x=x 17 ' 

n I n 
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then applying the boundary conditions, we obtain the recursion 
f ormulae t 


z in = z n 

n n 


\* Z nV anh q n d nJ 


(8) 


a'nd 


„ An + 1 + Cl\e- ( %- q n + l )x n (9) 


W = A n Z~ ^ + z^ 


where d n is the layer thickness, d n = x^x^. Finally, 
within each layer, the coefficients in Eq (4), B R and A fl are 
related by, 


(1-0) 


1G Equations (4-10) completely specify the temperature at every point in 
the material and will be used later in our calculations of the 
thermoelastic response. Before proceeding with that analysis, we 
should note that the gamier ophone experiment can also be analyzed 
with these results. That experiment is essentially a measurement of 
15 the surface temperature, T 0 = T^), which only depends on 
Q o and Z* n , i.e., 

Q c (U) 
o -in 
. 1 

. and Z[ n in turn depends ot» all of the other Z^ n ' s through the 
recursion relation (8). 
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Having a complete knowledge of the temperature, wc can now obtain the 
thcrmoelast i c response by solving the clastic wave equation, 

where <f> is the clastic displacement potential and y is the 
thermoelast ic constant defined below. The displacement, u, and 
stress, a , are obtained from ^ as follows, 

u = 

dx 


°xx = -P"V . (13) 


and 


a = a 

yy . zz 


where E is Young's modulus, v is Poisson's ratio, and a is the linear 
thermal expansion coefficient. Also in Eq (12), k is the elastic 
wave vector , 

\; «[$ N As]* (»> 

and y is a thermoelast ic constant, 

y - ir]7« . (is) 

The solution of Eq. (12) can be expressed as the integral, 

£(x) = y^x 1 g(x,x') 7(x') T(x') (16) 

where g(x,x') is the Green's function describing the elastic response 
at x to a fi-function source at x f . That is, g is a solution of the 
equat i on, 

2 

--§ + k 2 g = j(x-n') (17) 
dx 
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and satisfies whatever boundary conditions are imposed on the 
system. Since the purpose of this paper is to show the effects of 
thermal properties on the thermoelost ic response, we w.lll assume 
constant elastic properties throughout the material. The problem is 
5 then a simple one-dimensional elastic half-space problem with the 
well -known solution, ^ 

g(x>x .) = -.l-jV^hc-x'U e -ik|x + x.iJ (18) 

The second term on the right hand side of Eq. (18), + e ~ ik l x+x '| % 
ensures that the boundary conditions at x=0 are satisfied: 
lO + e ~** c i x+x ' corresponding to the rigid boundary condition, .u(0)=0 t and 
_ e -ikjx+x , j correS ponding to the free boundary condition, ^ xx (°)~°* 

In all practical applications, the frequencies are low enough that 
the elastic wave vector will always be much smaller than any of the 
tWnnal wave vectors, i.e., k « q n Fur lhermore> the lhermal v . ave , 

15 being highly damped, becomes insignificant after but a few thermal 
wavelengths, so that in Eq. (16), kx'ocl throughout the region in 
which the temperature is significant. Finally, in this work we are 
interested in the thermoelast ic response for large x, well beyond the 
region of any temperature change, that is, for x>x*. An excellent 
20 approximation to Eq. (18) is then, 

1 e , rigid boundary 

ik 

g(x,x") = { (19) 
- x f e~ ikx , free boundary 
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Since x'<*. 1/k, the rigid boundary condition gives rise to o n*uch 
larger thcrmoclast ic response then does the free boundary condition. 
This was previously noted by Whitc^ 5 , The physical Interpretation 
of this effect is as follows. The temperature at each point in the 
material generates two waves, g(x,x') f propagating to the right and 
left of the source. When the left-going wave strikes the boundary, 
x=0, it is reflected and interferes constructively or destructively, 
with the original right-going wave, depending on the boundary 
condition at x=0. For the rigid boundary there is no phase change 
upon reflection and since kx / «l, the two right-going waves are 
essentially in phase interfering in a constructive manner. On the 
other hand, the free boundary results in a phase change upon 
reflection and therefore destructive interference occurs. However, 
the interference is not completely destructive because of the small 
additional phase lag, 2kx' f between the reflected and unreflected 
waves. As we shall see, it is the small differential signal 
resulting from this phase-lag that produces the important 
depth-capability of thermoelast ic thermal-wave imaging. 

Using Eq. (19) for the Green's function in ?q. (16), we obtain for 
the thermoel as t ic response, 



t rigid 


free 


(20) 
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For a system of N thermal layers we can write 

N 

* = X>n (2D 
n=l 

where ^ is the contribution from the nth layer, 

- i kx 

" l£ e J x dx ' V x ' } • ri B id 

n-1 


(22) 


-" 7 e 
n 


-ikx 


J "dx' x* T n (x«) , free , 


and T n is oblained using Eqs. (4-10). From #(x) we can obtain 
either the displacement u or the stresses a xx * *yy» or c 
using Eq- (13) . 

-O 

11 Two- layer Exam ple 

The results of the Theory Section can best be appreciated with the 
consideration of a specific example. We consider the two-layer 
system shown in Fig. l r with the second layer much thicker than the 
^ first layer, and extending well beyond the range of any thermal 
waves „ 

Using Eqs. (8) and (11) we find the surface temperature T given by 


- ?° [ Z l * Z 2 tanh q l d 1 
Z l \ Z 2 + Z l tanh ""^i 5 " 


= I -± _± 1L I (23) 


From this, then, we can obtain the response of a gas-microphone 

20 system for the two-layer sample. Note that Eq. (23) is exactly 

equivalent to the expression for T Q obtained from the 

27 

Rosencwa ig-Gersho theory for the type of sample considered here. 
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Next, from Eqs. (4-10) and (21-22), we obtain the thermoelastic 
response. 


for the case of a rigid front surface boundary, and 


(24) 


7Q„ -ikx 


w/>C 


"»1 



1 + tanh q^d 


1 + li tanh q,d 


(25) 


for the case of a free front surface boundary. Equation (24) is 

15 

equivalent to White's result for a rigid boundary. 


Note that we are here considering a sample with an inhomogeneity only 
in the thermal conductivity . For such a sample, Eq. (24) shows 

lO that the thermoelastic response -when the front surface is rigidly 

clamped will be independent of at , and thus no information about the 
second layer can be obtained for such a sample under rigid boundary 
conditions. When the front surface (heated surface) is free, the 
thermoelastic signal will now be much smaller, but it will contain 

15 information about the second layer. Thus, only the surface 

temperature T o and the free surface thermoelastic signal^ provide 
information about the second layer. 


6AD ORIGINAL 


0097473 


- 36 - 


We should point out at this time that Wetsel, in his treatment , 
has developed a similar two-medium one dimensional model for 
thermoelast ic wave generation where the emphasis is cn the efficiency 
of elastic wave generation rather than on imaging of subsurface 
5 thermal features. In his work, Wetcel has shown that air and vacuum 
at the front interface of the sample produce quite different 
results- In the analysis presented here, we assume a vacuum 
condition only. A future paper will discuss our analysis with 
respect to air at the interface. 

lo Let us now consider two limiting cases, one at low frequencies where 
v:e have long thermal wavelengths and d/p < 1, and the other at high 
frequencies where we havp short thermal wavelengths and 6/p > 1. 
-For the long thermal wavelength limit with d//i < 1, we find 
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In .the Ihcrmoelast ic response we find 


<A ( x ) = i_e S. 


( qi d) 2 /zJ 


(27) 


for the free boundary case* Equation (24) holds for the rigid 
boundary case. 


L5 


As we said before, the rigid boundary thermoelas t ic signal provides 
no information about the second layer, and only the surface 
temperature T o and the free surface thcrmocl as t i c signal 0 provide 
such information. Note also that in the long thermal wavelength 
case, T e varies linearily with d/p while <p Q varies quadrat ically 
with d/// . In general the information that will be desired are the 
.distance d between the two layers and the thermal conductivity * 9 
of the second layer (or the ratio Z^/Zj) . In examining Eqs. (26) 

.and (27) we see that unambiguous and unique evaluations of d and * 2 
would be very difficult in the long thermal wavelength limit since 
here both the magnitude and phase of the signal^ are determined by 
terms that are products of d nnd 


20 


In the short thermal wavelength limit d/j/ > 1 we find 
and 


_7Q 0 e -ikx 

' upC ~Tq~ 


1 + 2 s 


-.qua t ion (24) again applies for the rigid boundary case, 


(28) 


(29) 
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Unlike the long thermal wovelcnglh cose, the shor t - thermal wavelength 
situation provides on oppor tun i ty for an unambiguous determination of 
both d and * 2 by measuring magnitude and phase. This is possible 
since there are now ho phase terms that are products of both x 0 and 
d, and in fact, there are no phase terms that are functions of 
* 2 . Thus, the real parts of Eqs. (28) and (29) can be written as 

= T; £l + JlU 2 )e- 2d//, l cos( ^ - ad/y/^j (30) 


T = 
o 


and 


*,(x) = ^J(x) £l + TU 2 )e d/ ^l cost 0, - d/y/^J (31) 

In these equations T^ and represent the signal amptitudes that 

we would obtain with no second layer present, and ^ and ip 2 the 

corresponding phases;, The quantity H( * 2 )e~ 2d,;, l represents the 

"correction to the T^ term arising from the reflection of the 

thermal wave from the second layer, while the K^Je^^l term 

represents the correction to the^(x) term arising from the 

transmission of the thermal wave into the second layer. These 

magnitude corrections are functions of both * 2 and d. On the other 

hand, the phase correction terms are functions of d only. Thus, by 

measuring magnitude and phase both k % and d can be determined 

uniquely. If one has N layers in the sample, then N frequencies 

would be needed to define the entire set of d and k 

n n" 

In Figs. 2 and 5, wc present some graphical illustrations of the 
surface temperature and thermoclas t ic signals by showing how the 
total magnitude and phase of the T o and £ ff signals* vary with 
frequency for the two cases of * 2 /lf i = 4 and *2 /Ar i = °- 25 - 
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The signals have all been normalized to a single layer signal, that 
is, we plot M 2 / M i where M 2 and are the magnitudes of 
either the surface temperature or the Ihermoel as t ic response for the 
two-layer (Mg) and single-layer (M^ cases* The symbol 5 T Q 

5 and $ 0 identify the curves corresponding to the front surface 

temperature and free surface thermoel as t i c signals respectively. The 
phase has also been normalized by subtracting the phase for the 
single-layer case, i.e., \p= ^ - The frequency range is also 

normalized and runs from d//^ = 0.1 to d//^ = 10. We can see 

10 from Figs. 2 and 5, that the thermoelastic response is more sensitive 
to the presence of the second layer than the surface temperature 
response. In addition, the presence of thermal-wave interference as 
evidenced-by the oscillations in the magnitudes, when d/j/j is 
between 1 and 3, is much more evident in the thermoelastic signal. 
15 Both of these effects crise from the fact that the thermoelastic 
signal from the seeond layer is attenuated by e'^^l while for the 
surface temperature, the second layer signal undergoes a much greater 
e 2d ^i attenuation. This is, of course, a result of the 
transmission nature of the thermoelastic response compared to the 
20 reflection nature of the surface temperature response. 

\ * v * Subsurface Heating 

.J Our treatment to this point has assumed that ell of the energy 
i absorption and initial heating occurs at the front surface* This is, 

-25 ot course, somewhat unrealistic since energy absorption and initial 
heating r/ill always occur down to some depth beneath the surface. 
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The case of subsurface heating is an important one, not only because 
it is more realistic, but also because it has important implications 


for t hcnnal-wave resolution and image depth. For example, in some 

recent work on thermal -wave electron imaging, Cargill has 
14 

suggested- that the image resolution is set not only by the thermal 
wavelength and the size of the energy beam at the surface, but also 
by the total beam broadening that occurs as the electron beam 
penetrates into the sample* 

To examine this situation we consider a simple three-layer sample as 
shown in Fig* 6 where the heating is assumed to occur at a subsurface 
layer x=x q , and that the distance between the front surface, x=0, . 
and x=x c is_ d . We already know that only a stress-free front 
surfece will produce thcrmoelast ic signals that provide thermal -wave 
images due to variations in thermal conductivity. Thus we will 

-assume in the sample of Fig. 6 that the front surface is 
stress-free . 

Proceeding as before in the case of heating at the surface, we first 
find the temperature at every point within the material. For 
subsurface heating, T(x) is still given by Eq.(4) with Eqs.(8-10) 
valid for n>l . For n=0 v/e have 



cosh o d ' e H l o 


o o 


(32) 


with 


A 


Q c sech q Q d 0 


(33) 
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for the surface temperature, T*, we obtain 


T S = T e'Vo 
o o 


1 + tonh q„d rt 
o o 


1 + r-~ lanh q d 
„in < 

Z l 


(35) 


where T Q is the result obtained earlier for heating at the surface 

given by Eq. (11). This result of course reduces ip T when 

o 

d o =0, and goes to zero for large d Q . 

Knowing the temperature, we can now proceed exactly as before using 
the Green's function given by Eq. (19) in the expressions for the 
therxnoelast ic. response, Eqs. (20-22). Applying this analysis to the 
three-layer example shown in Pig. 6, we obtain 


1 4 -7 tanh q d 
_in o o 


1 " v»,V^ 



1 - sech q d 
_o_o_ 

Z 

1 + _°.tanh q d 
«in © o 


(36) 

where £§(x) denotes the thermoelas t ic response for subsurface 
heating and^(x) is the result for the two-layer surface heating 
problem in Fig- 2 given by Eq.(25).The second term in Eq- (36) is a 
small" correction term as long as d Q < p Q However when 
d 0 >/i 0> the second term becomes comparoble tq and then quickly 

-1 

larger than the first term. Finally when d Q » p but still <k 
t the thermoelas t ic signal is given by 
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This expression is essentially identical to that for the rigid 
boundary case, Eq. (24), except that d Q is used in place ofl/k • 
Thus 0^(x) will be smaller in magnitude than 0 u (x), but just as 
5 with 0 o (x), #^(x) will also be completely independent of the thermal 
conductivity. 

The results summarized inEqs.(36) and (37) arc most significant. 
They indicate that although energy absorption and heating beneath the 
surface will produce thermoelastic signals, these signals will not 
D ■ image any variations in thermal conductivity when the heating occurs 
more than a few thermal diffusion lengths beneath a stress-free 
mechanical surface. This in turn implies that where thermal-wave 
Images are dominated by variations in the thermal conductivity, such 
images will have a resolution determined by the thermal wavelength 
15 j»nd by the size of the incident beam during its propagation within a 
theamal wevelcngth beneath the surface. The image depth is also 
limited to only a few thermal diffusion lengths (—one thermal 
wavelength) beneath the surface. This is a result independent of the 
total depth of penetration of the heating beam, and thus the 
20 resolution is independent of the amount of beam spreading that occurs 
at depths beyond the first few thermal diffusion lengths. 
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It is important to note that the only assumption made in deriving 
this result is the reasonably realistic one that the major thermal 
variation in the sample is in its thermal conductivity. Note also 
that in the limit ^ o =0 > we recover our earlier results for heating 
5 occurring at the front surface. 

V. Conclusi ons 

We have derived a model for thermal-wave depth-profiling that is 
applicable both to conventional gas-microphone and photothermal 
thermal-wave measurements where the measured parameter is the 
lO periodic temperature at £he surface of the sample, and to 

piezoelectric Rhotoacoust ic and thermal-wave experiments where the 
signal is the result of a su^suirface thermoe) us t ic response. We have 
shown that the significant differences in thermal-wave imaging 
capability between samples with mechanically free and constrained 
15 surfaces, first predicted'*by Whi te 15 , for thermoclast ic ultrasonic 
generation only, can be analyzed in terms of wave interference 
effects. Mechanically constrained surfaces produce greater 
thermo-elast ic signals but these signals are independent of 
variations in thermal conductivity, and thus cannot generate 
20 thermal-wave images of such variations. Mechanically free surfaces 
produce smaller signals but these signals do provide images of 
variations in thermal conductivity. V/e also find that the 
thermoelast ic signal is a more sensitive measure of any subsurface 


* 
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variations in thermal conductivity than the surface temperature 
signal* 

We have used as a model a sample that has nonuniform thermal but 
uniform elastic properties since it provides a simpler analysis and a 
cleaner interpretation." It also clearly shows both the inherent 
similarities and the differences between thermal-wave imaging 
performed by measuring surface temperatures, as in gas-microphone and 
photothermal experiments, and thermal-wave imaging performed by 
measuring the thermoelast ic response in the sample. A model of a 
■ LO sample with nonuniform thermal and uniform elastic properties is* 
moreover appropriate for the most important application of 
thermal^ave depth-profiling of dopants in semiconductors. 
Por such a model at any given modulation frequency, only two unknown 
parameters, such as a length and a thermal conductivity, can be 
15 derived from amplitude and phase measurements of either the 
thermoelast ic or thfe surface temperature signals. When more 
parameters need to be determined, more frequencies must be used, and 
if an accurate depth profile is to be obtained thermal-wave signals 
at several frequencies must be measured. Furthermore, unambigious 
20 determinations of layer thicknesses and thermal conductivities can 
best be obtained at frequencies that produce thermal diffusion 
lengths, shorter than the thicknesses to be determined. 

We have also considered the case where the heating may occur beneath 
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the surface of the sample. We f\nd that z»s long as the heating 
occurs within one thermal diffusion length beneath the sample 
surface, the results are not significantly different from the case of 
front surface heating. However, for heating occurring at greater 

5 depths, the thermoelastic signals become independent of variations in 
thermal conductivity and thus will not provide thermal-wave images of 
such variations. Thus in thermal-wave electron or laser microscopy, 
the imaging depth is within a few thermal diffusion lengths of the 
surface even though the electron or laser beam may penetrate much 

lO deeper. Furthermore, thermal-wave resolution will not be degraded by 
any beam spreading that occurs beyond this thermal-wave imaging depth 

Finally, we point out that althpugh the model in this paper is 

one-dimensional, we expect that, just as in photoacoustic theory , 

three-dimensional treatments will not significantly alter the 

15 physical understanding of thermal-wave imaging derived from this 

one-dimensional treatment. It is true that edge effects can produce 

signal changes and that these edge effects can only be analyzed with 

a three-dimensional theory. However, an analysis of these edge 

effects indicates that their affect on the signal is a function of 

20 the ratio of the thermal diffusion length to the beam spot size. 

o q 

Thus in the photothermal work of Murphy and Aamodt , the edge 
effects are particularly strong because, at the low modulation 
frequencies used, the thermal diffusion length is much greater than 
the diameter of the focused beam spot. On the other hand, at the 
25 high modulation frequencies employed In the high-resolution 
thermal-wave microscopy of interest here, the thermai diffusion 
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lengths are comparable to or smaller than the beam spot size and the edge 
effects can be expected to be relatively minor. This appears to be true 
in fact, since in none of the high-resolution experiments reported by 
Rosencwaig and his colleagues 12 * 13 * 22, 30 " 32 , or by Cargill 14 * 33 , are 
there any indications of noticable edge effects. Thus, we believe that 
this one-dimensional theory provides an adequate, if not complete, 
representation of the physical processes occurring during thermal-wave 
imaging. 
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CLAIMS 

1* A method for evaluating the thickness of a layer 
or layers of material on a substrate by measuring 
5 either the phase or magnitude parameters of thermal 
wave signals generated by a focused periodic heat 
source, characterised by the steps of: 

focusing said periodic heat source on the 
uppermost layer deposited on said substrate; 
10 measuring the value of one of said parameters 

of the thermal wave signals generated in said 
sample at at least one selected modulation 
frequency of said source; and 

comparing the obtained value of said one 
15 parameter with a predetermined value of said one 

parameter associated with a reference sample 
whereby the thickness of the layer can be 
evaluated. 

20 2. A method according to claim 1, characterised in 
that said substrate is a semiconductive material such 
as silicon* 

3. A method according to claim 1 # characterised in 

25 that the thermal wave signals are measured with a means 
that detects the oscillating temperature of the heated 
spot on the surface of the substrate. 

4. A method according to claim 3, characterised in 

30 that said detection means includes a photoacoustic cell 
apparatus. 

5. A method according to claim 3, characterised in 
that said detection means is defined by a laser Which 

35 senses the periodic heating of a medium in contact with 
the heated spot on the surface of the substrate. 
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6. A method according to claim 3, characterised in 
that said detection" means is defined by an infrared 
3ete£tor which senses the periodic infrared emissions 

from the heated spot on the surface of the substrate. 

mxG'ttit y; • . - i 1 

7. , A method according to claim 1, characterised in 
that said thermal wave signals are measured with a 

, 4 means that detects the oscillating thermal displacement 
of the surface of the substrate at the heated spot* 


IP 


8. A method according to claim 7, characterised in 
that, js^Ld detection means includes a laser probe* 

t- - :: - # - 

J9_._ ^^jn^thod accprdit^g to claim 7, characterised in 

15. that jsaid detection means includes a laser 

interferometer . 

10. A method according to claim 1, characterised in 
that the thermal wave signals are measured with a means 
20^. .for detecting the acoustic signals that are generated 
by the thermal waves. 

■ 11. A method according to claim 10, characterised in 
that said detection means includes an ultrasonic 
25 transducer. 

12. A method according to claim 10, characterised in 
that said detection means includes a laser probe. 

30 13. A method according to claim 10, characterised in 
that said detection means includes a laser 
interferometer - 

14. A method for determining the thickness of a layer 
35 or layers of material deposited on a substrate by 

measuring either the I phase or magnitude parameters of 
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thermal wave signals generated by a focused periodic 
heat source, characterised by the steps of: 

focusing said periodic heat source on the 
uppermost layer deposited on said substrate; 
5 measuring the value of one of said parameters 

of the thermal wave signals at a plurality of 
selected modulation frequencies of said source 
wherein the number of modulation frequencies 
selected is greater than the number of layers 
10 whose thickness is to be determined; 

normalizing the measured value of said 
parameter of each modulation frequency selected 
relative to the value of said parameter determined 
for a reference sample; and 
15 comparing said normalized values with expected 

normalized values derived from a model depicting 
the thermal process in said samples vfoereby the 
thickness of said layers can be determined. 

20 15. A method according to claim 14, characterised in 

that the thermal wave signals are measured with a means 
that detects the oscillating temperature of the heated 
spot on the surface of the uppermost layer deposited on 
the substrate. 


25 


30 


16. A method according to claim 15, characterised in 
that said detection means includes a photoacoustic cell 
apparatus • 

17. A method according to claim 15, characterised in 
that said detection means is defined by a laser which 
senses the periodic heating of a medium in contact with 
the heated spot on the surface of the uppermost layer 
deposited on the substrate. 


35 
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18. A method according to claim 15, characterised in 
that said detection means is defined by an infrared 
detector which senses the periodic infrared emissions 
from the heated spot on the surface of the uppermost 


19. A method according to claim 14, characterised in 
that said thermal wave signals are measured with a 
means that detects the oscillating thermal displacement 

10 of the surface of the uppermost layer deposited on the 
substrate at the heated spot. 

20. A method according to claim 19, characterised in 
that said detection means includes a laser probe. 


21. A method according to claim 19, characterised in 
that said detection means includes a laser 
interferometer • 

70 22. A method according to claim 14, characterised in 

.that the thermal wave signals are measured with a means 
which detects the acoustic signals that are generated 
by the thermal waves. 

25 23. A method according to claim 22, characterised in 
that said detection means includes an ultrasonic 
transducer. 

24. A method according to claim 22, characterised in 
30 that said detection means includes a laser probe. 

25. A method according to claim 22, characterised in 
that said detection means includes a laser 
interferometer . 


5 


layer deposited on the substrate. 


15 


35 
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26. A method according to claim 14, characterised in 
that the parameter measured is the phase of the thermal 
wave signals and the values are normalized by 
subtracting the value of the phase of the reference 
sample from the phase of the multilayer substrate. 

27. A method according to claim 14 , characterised in 
that the parameter measured is the magnitude of the 
thermal wave signals and wherein the values are 
normalized by taking a ratio of the value obtained from 
the reference sample and the value obtained from the 
multilayer substrate • 

28. A method according to claim 14, characterised in 
that the comparison of the measured normalized values 
to the normalized values derived from the model is done 
by a least-squares fitting routine. 

29. A method for determining the dtiange in thermal 
characteristics as a function of depth in a non-uniform 
sample having impurities or defects therein by 
measuring either the phase or magnitude parameter of 

"thermal wave signals generated in said sample by a 
focused periodic heat source, characterised by the 
steps of; 

focusing said periodic heat source on said 
non-uniform sample; 

measuring the value of one of said parameters 
of said thermal wave signal generated in said non- 
uniform sample at a plurality of seleted 
modulation frequencies; 

normalizing the measured value of said 
parameter at each modulation frequency relative to 
the value of said paramter determined for a 
reference sample; and 

comparing said normalized values with 
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expected normalized values derived from a model 
depicting the thermal process in said samples , 
with said model for the non-uniform sample being 
characterized as divided into a plurality of 

5 hypothetical layers each having the same 

thickness, with the number of layers being less 
than the number of modulation frequencies 
selected, whereby the thermal characteristics of 
the hypothetical layers and the depth variation of 

10 these thermal characteristics in the non-uniform 

sample can be determined. 


30. A method according to claim 29, characterised in 
that the thermal characteristics are the thermal 

15 conductivities of the hypothetical layers and further 
including the step of correlating the concentration of 
impurities or defects in the sample to variations in 
the thermal conductivity such that a profile of the 
impurities or defects can be obtained as a function of 

20 depth. 

31. A method according to claim 29, characterised in 
that the measurements of said one parameter are taken 
at a plurality of modulation frequencies whereby the 

25 resolution of said depth variations is increased. 

32. A method according to claim 29, characterised in 
that the thermal wave signals are measured with a means 
that detects the oscillating temperature of the heated 

30 spot on the surface of said non-uniform sample. 

33. A method according to claim 32, characterised in 
that said detection means includes a photoacoustic cell 
apparatus . 
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34. A method according to claim 32, characterised in 
that said detection means is defined by a laser which 
senses the periodic heating of a medium in contact with 
the heated spot on the surface of said non-uniform 

5 sample* 

35. A method according to claim 32, characterised in 
that said detection means is defined by an infrared 
detector which senses the periodic infrared emissions 

10 from the heated spot on the surface of said non-uniform 
sample. 

36. A method according to claim 29, characterised in 
that said thermal wave signals are measured with a 

15 means that detects the oscillating thermal displacement 
of the surface of said non-uniform sample at the heated 
spot. 

37. A method according to claim 36, characterised in 
?0 that said detection means includes a laser probe. 

38. A method according to claim 36, characterised in 
that said detection means includes a laser 
interferometer . 

25 

39. A method according to claim 29, characterised in 
that the thermal wave signals are measured with a means 
which detects the acoustic signals that are generated 
by the thermal waves. 

30 

40. A method according to claim 39, characterised in 
that said detection means includes an ultrasonic 
transducer. 

35 41. A method according to claim 39, characterised in 
that said detection means includes a laser probe. 
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42- A method according to claim 39, characterised in 
that said detection means includes a laser 
interferometer . 

43. A method according to claim 29, characterised in 
that the parameter measured is the phase of the thermal 
wave signals and the values are normalized fcy 
subtracting the value of the phase of the reference 
sample from the phase of the non-uniform sample. 

44. A method according to claim 29, characterised in 
that the parameter measured is the magnitude of the 
thermal wave signals and wherein the values are 
normalized by taking a ratio of the value obtained from 
the reference sample and the value obtained from the 
non-uniformed sample. 

45. A method according to claim 29, characterised in 
that the comparison of the measured normalized values 
to the normalized values derived from the model is done 
by a least squares fitting routine. 
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